Low energy transition in spectral statistics of 2D interacting fermions 
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to investigate efficiently the transition from nonergodic 
(localized) to ergodic eigenstates. 

Here we use the above method to study the change of 
the spectral statistics, P(s), with excitation energy E in 
a model of spinless fermions with Coulomb interaction 
on 2D disordered lattice. The Hamiltonian reads 
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We study the level spacing statistics P(s) and eigenstate properties of spinless fermions with 
Coulomb interaction on a two dimensional lattice at constant filling factor and various disorder 
strength. In the limit of large lattice size, P(s) undergoes a transition from the Poisson to the 
Wigner-Dyson distribution at a critical total energy independent of the number of fermions. This 
implies the emergence of quantum ergodicity induced by interaction and derealization in the Hilbert 
space at zero temperature. 
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The experimental observation of the metal-insulator 
transition in two dimensions (2D) by Kravchenko et al. 
has attracted a great interest to interacting fermions 
in a disordered potential. Indeed, according to the well- 
established result 0], all states are localized for nonin- 
teracting particles in 2D. Therefore, in the view of the 
experimental result [jjj, a new theory should be devel- 
oped to understand the interaction effects between the 
localized fermionic states. However, in spite of various 
theoretical attempts, a coherent theory for such systems 
is still not available. While for highly excited states, it 
has been shown that the repulsive/ attractive interaction 
can induce a derealization of two interacting particles 
[^-^| , the properties of low energy states are not under- 
stood yet. Recently, in addition to experimental and the- 
oretical investigations, a number of attempts have been 
made to study these many fermionic systems through 
numerical simulations |7|-[|. Even though several inter- 
esting features have been reported, the systems studied 
there are not large enough to observe interaction induced 
derealization. 

In this paper we use another numerical approach based 
on the analysis of spectral properties of multi-particle 
fermionic systems. Indeed, Shklovskii et al. have shown 
that the level spacing statistics is a powerful tool to an- 
alyze the Anderson transition in disordered systems [flO| . 
When the states are localized, the levels are not corre- 
lated and the statistics is given by the Poisson distri- 
bution, P(s) = Pp(s) = exp(— s), while in the metallic 
phase, the states are ergodic and the statistics is close to 
the Wigner surmise, Pw(s) = (its/ 2) exp(— tts 2 /A). The 
critical transition point is characterized by an intermedi- 
ate statistics which depends on the boundary conditions 
p[ and the spatial dimension of the system [Q . This 
approach has also been used to determine the quantum 
chaos border and the interaction induced thermalization 
in finite fermionic systems |l3[ ] and to detect Anderson 
transition for two electrons with the Coulomb interac- 
tion on 2D disordered lattice JT^ , p^[ . All these results 
demonstrate that the approach developed in pQ] allows 



where aj( a ») is the fermion creation (annihilation) op- 
erator at site i, V is the hopping between the nearest 
neighbors, the diagonal energies Wi are randomly dis- 
tributed within the interval [—W/2,W/2], m = a\a,i is 
the occupation number at site i and U is the strength 
of the Coulomb interaction with nj the interparticle dis- 
tance. The particles move in a 2D cell of size L x L with 
the periodic boundary conditions applied. The Coulomb 
interaction is taken between electrons in one cell of size 
L and with 8 charge images in nearby 8 cells as in juj. 
The number of particles N p and the cell size were varied 
within the intervals 2 < N p < 20 and 8 < L < 28. 
With the notations in the Eq. (1), the parameter r s , 
which measures the ratio of the Coulomb energy to the 
Fermi energy, ep, is given by r s = U / (2V \pKv) , where 
v = Np/L 2 is the filling factor and tp = iiri/V. The 
majority of our data have been obtained for v w 1/32 
(nearest rational value) and for U/V = 2, which corre- 
sponds to r s = 3.22. 

To study the level spacing statistics P(s), we general- 
ize the approach used in |l4f| for many particles. First, 
one particle eigenstates (orbitals) at U = are obtained 
and the Hamiltonian (1) is rewritten in this basis us- 
ing the two-body interaction matrix elements. We con- 
sider the first M orbitals from the lowest energy and 
the Hamiltonian multi-particle matrix, constructed only 
from these orbitals, is diagonalized at the final stage. To 
the maximum, N p = 20 particles with up to M = 42 or- 
bitals has been considered and the resulting matrix size is 
N m 5000. This size is significantly smaller than mCn 



since the condition ^ 



mi < Y^i=i 1 * + M applies for 
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the one-particle orbital index m,. Such a pyramid rule 
allows one to use efficiently only low energy multi-particle 
states and to make a striking reduction of the resulting 
total matrix size without any serious modification of the 
low energy states. We checked that our level statistics 
results at low energy are not sensitive to the variation of 
M and N m (see for example insert in Fig.l). To compute 
the spectral statistics P(s) at a given total excitation en- 
ergy E, which is counted from the ground state energy, 
disorder average has been performed over Np> = 5000 
(for low energy) and Np, — 1000 (for higher energy) con- 
figurations. In this way, the total statistics for P(s) for 
small energy interval varies from N$ = 10 4 at low E to 
N s = 3 x 10 5 at high E. 




FIG. 1. Level statistics P(s) for 10 particles at L = 18, 
r s = 3.22 in the energy interval 0.25 < E/B < 0.3 (with 
B = 4V): W/V = 15,r) = 0.93 (o); W/V = 10, rj = 0.51 (o); 
W/V = 7, rj = 0.10 (•). Total statistics is iV s = 10 5 ; 1.4xl0 5 ; 
1.8 x 10 5 , respectively. Full lines show the Poisson distri- 
bution and the Wigner surmise. Insert, made for N p = 6 
case of Fig. 2b, illustrates that at low E/B the statistics P(s), 
characterized by rj, is independent of matrix size variation 
N m = 797(o), 1231 (□), 2235 (full diamond). 



We note that for given values of disorder the one-particle 
inverse participation ratio (IPR), £i (the number of sites 
contributing a one-particle eigenstate) is much smaller 
than the total number of sites L 2 for the case of Fig. 1: 
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A change of P(s) with W, at a given E and all other 
parameters fixed, is shown in Fig. 1 for N p = 10. As 
W decreases, P(s) evolves from the Poissonian to the 
Wigner-Dyson distribution. To measure the proximity of 
P(s) to one of these two limiting cases, it is convenient to 
use the parameter r\ which is defined as r\ = J Q ° (P(s) — 
P w {s))ds/ Jo°(Pp(s) - P w (s))ds, where s = 0.4729... 
is the smaller intersection point of Pp(s) and Pw(s). In 
this way 77 = 1 corresponds to Pp(s) and 77=0 to Pw(s). 
According to Fig. 1, rj changes almost by an order of 
magnitude when W decreases by factor of two. This 
shows that for strong disorder, W/V — 15, the multi- 
particle states at given E are not ergodic (localized) while 
for weaker disorder, W/V = 7, they become ergodic and 
characterized by the random matrix spectral statistics. 



FIG. 2. Dependence of 77 on the rescaled total energy E/B 
for various number of particles N p = 2(A), 3 (full triangle 
up), 4 (□), 6 (full diamond), 8 (o), 10 (•), 14 (o) and 20 
(full triangle down); (a) W/V = 10 and (b) W/V = 7; filling 
factor v ~ 1/32 and r s — 3.22, 8 < L < 25. Inserts show the 
dependence of £s/N p on E/B by same symbols; arrows mark 
the critical energy E c from the main figure. 

6 = 3.4 and 4.2 (W/V = 15); 5.2 and 11.6 (W/V = 
10); 8.2 and 36.7 (W/V = 7), where the formers are for 
the ground state and the latters for the center of the 
band. This means that the transition to ergodicity is in- 
duced by the interaction which becomes effectively more 
strong when the one-particle localization length increases 
(similar to the two interacting particle case |-§). At the 
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same time, the comparison with the data for N p — 2 H] 
at the same E (r) « 0.92, 0.82 and 0.51 for W/V = 15, 
10 and 7, respectively) shows that the derealization ef- 
fect due to multi-particle interaction is stronger for weak 
disorder while for strong disorder it does not affect local- 
ization. 



P(s) 




FIG. 3. Level statistics P(s) at the crossing points of Fig. 2 
(with the same symbols). The upper data are for W/V = 10, 
0.225 < E/B < 0.275, rj c w 0.56, iV s = 10 4 ,4 x 10 4 and 10 5 
for N p = 3, 6 and 10, respectively. The lower data (shifted 
down by 0.4) are for W/V = 7, 0.125 < E/B < 0.175, 
r] c x 0.33, N s = 10 4 ,1.4 x 10 4 and 5 x 10 4 for N p = 3, 6 
and 14, respectively. The full lines are the Poisson distribu- 
tion and the Wigner surmise. 

In Fig. 2 we present the variation of r\ with the to- 
tal energy E for various number of particles at fixed 
filling factor and two values of disorder. For N p > 2, 
all curves have approximately the same crossing point: 
E c w 0.255 (W/V = 10), E c w 0.15B (W/V = 7) and 
E c w 0.1B (r] c w 0.19, W/V = 5, not shown). As N p 
increases, the statistics approaches the Poissonian case 
for E < E c , while it goes to the Wigner case for E > E c . 
The transition point, E c , is characterized by an interme- 
diate statistics, which is independent of the system size 
as shown in Fig. 3. The energy E c grows with the in- 
crease of W so that for W/V = 15 the crossing point is 
not found clearly within our computations (in this case 
X) fa 0.8 - 1.0 for < E/B < 1 and N p > 2). The above 
data give a strong evidence that the ground state and 
the excited states with E < E c are localized in the limit 
of the infinite system size. However, for excitations with 
E > E c the statistics is close to Pw(s) that implies the 
ergodicity of eigenstates and their space delocalization. 

To show that the eigenstate properties are qualitatively 
different below and above E c we determined the varia- 
tion of IPR £s with N p and E as presented in the inserts 



of Fig. 2. This IPR £s is computed in the basis of nonin- 
teracting Slater determinants and determines how many 
of such states are required to construct an eigenstate. 
The results are obtained by averaging over 100 disorder 
realizations so that the statistical error is less than 5%. 
For E < E c the ratio £s/-/V p , which gives the number of 
nonintcracting states per particle, remains constant while 
for E > E c it grows rapidly with N p . For example, the 
ground state average value is £s/N p = 0.35; 0.40; 0.46 
for W/V = 10; 7; 5 (the latter is not shown) with only 
±21; 13; 18% variation when N p changes from 3 to 10; 
14; 20, respectively for each W/V. On the contrary, for 
E = 0.45 > E c the ratio £s/N p grows in 2.7; 9; 20 times 
when N p changes from 3 to 8 for W/V = 10; 7; 5 re- 
spectively. These results show that the delocalization in 
the Hilbert space of Slater determinants takes place for 
E > E c with E c independent of number of particles. A 
crossover from localized to extended states in the Hilbert 
space was extensively discussed recently for the metallic 
phase with extended one-particle states [^,^3) . Our data 
in Fig. 2 represent the first evidence for a delocalization 
transition in the Hilbert space for the insulating phase 
with localized one-particle states. 
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FIG. 4. Dependence of e v /B on the number of particles 
N p , obtained from Fig. 2: W/V = 10 with r) (E v ) = 0.4 
(full diamond) and W/V = 7 with 77 (E v ) = 0.2 (•), 
where e v = E v /N p . The straight line shows the slope when 
E v — const. Insert gives the dependence of maximal r\ on 
r s for W/V = 7 and N p = 6: U/V = 2, 8 < L < 28 (full 
diamond) and L = 14, 0.25 < U/V < 2 (O). 

An unusual property of the above transition is that it 
takes place at the finite total energy. This means that 
the energy per particle e = E/N v , or the temperature, is 
equal to zero at the transition point as N p — > 00. In this 
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sense this transition can be considered as a quantum zero 
temperature transition. To ensure that in our numerical 
computations e was sufficiently low, we present variation 
of e v , defined at the fixed r\ level, with the number of 
particles in Fig. 4. According to these data e n drops 
more than by an order of magnitude and becomes ap- 
proximately by an order of magnitude smaller than the 
Fermi energy tp ps 0.1 B. 

Dependence of rj m on r s is found in the insert of Fig. 4. 
The value of r s is varied in two different ways: i) by 
changing the cell size L or ii) by changing U/V. The 
first case i) allows to obtain large r s value and shows that 
the statistics approaches the Poissonian case. This result 
is in a qualitative agreement with the experimental data 
0, where the localized (insulating) phase appears for 
large r s , and with the theoretical argument given in [[l4"[ , 
according to which the two-body interaction drops as the 
filling factor decreases. The same arguments explain why 
ij goes to 1 when the interaction strength U decreases as 
in the second case ii). One can expect that the variation 
of rj c with r s at large N p and fixed v will be qualitatively 
similar to that one of rj m in Fig. 4. 

The comparison of the energy E c , at which the transi- 
tion at fixed v takes place (see Fig. 2), with the transition 
energy for two particles (Ei c jB = 1.2 for W/V—10 and 
E 2c /B = 0.56 for W/V=7 @) shows that E c is signif- 
icantly smaller than E% c (E c jEi c ~ 0.2). This signifies 
that the multi-particle interaction is more efficient for de- 
localization than the two-body interaction for only two 
particles. Our qualitative understanding for the transi- 
tion to the quantum ergodicity at fixed v and fixed total 
energy E c is based on the following scenario. In the anal- 
ogy with the approach developed in pi] for N p — 2, the 
dynamics of three or four particles at finite energy E can 
be assumed to be equivalent to one particle dynamics 
in a disordered system with effective spatial dimension 
3 < d e ff < 2N p . With the growth of E, the rate of this 
one-particle hopping increases ju[ and the derealization 
transition finally takes place. As a result, in the original 
system, the energy can be transferred by this small group 
of particles from one place to another that allows to es- 
tablish the ergodicity in the system of large size at fixed 
filling factor. The values of r\ c found here (Fig. 2) are 
larger than that for 3D Anderson model [r\ c w 0.2), that 
makes the above scenario consistent with the result of 
fL2| according to which r\ c grows with the increase of the 
spatial dimension d. We also note that the approach of 
rj to 1 near the ground state is not so surprising. Indeed, 
even in the metallic quantum dots one has 77 rs 1 near 
the ground state since the interaction between particles 
is effectively weak JlJ,[lf| and the transition to ergodicity 
with 77 = takes place only at higher energy p7| . 

In conclusion, our results show that 2D fermions, which 
are strongly localized by disorder without interaction, be- 
come ergodic due to Coulomb interaction. The transition 
to ergodicity, in lattice and Hilbert spaces, and the ran- 



dom matrix statistics takes place at the finite total energy 
or zero temperature. These results are in favor of the ex- 
perimentally observed metal- insulator transition [Q , even 
if our approach does not allow to investigate the con- 
ductance dependence on temperature which is the main 
experimental method to detect the transition. Further- 
more, they show that the variable range hopping trans- 
port can be induced by electron-electron interaction in 
agreement with recent experiments [p| . 
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